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Abstract 

We explore the possible cross sections for the elastic scattering of neutralinos x on nucleons 
p,n in the minimal supersymmetric extension of the standard model (MSSM). Universality 
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of the soft supersymmetry-breaking scalar masses for the Higgs multiplets is not assumed, 
but the MSSM parameters are nevertheless required to lead consistently to an electroweak 
vacuum. We explore systematically the region of MSSM parameter space where LEP and 
other accelerator constraints are respected, and the relic neutralino density lies in the range 
0.1 < Q x h 2 < 0.3 preferred by cosmology. We also discuss models with Q x h 2 < 0.1, in which 
case we scale the density of supersymmetric dark matter in our galactic halo by Q x h 2 /0.1, 
allowing for the possible existence of some complementary form of cold dark matter. We 
find values of the cross sections that are considerably lower than the present experimental 
sensitivities. At low neutralino masses, m x < 100 GeV, the cross sections may be somewhat 
higher than in the constrained MSSM with universal soft Higgs masses, though they are 
generally lower. In the case of large m x , the cross sections we find may be considerably larger 
than in the constrained model, but still well below the present experimental sensitivity. 
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1 Introduction 



One of the key issues at the frontier between particle physics and cosmology is the nature of 
the non-baryonic dark matter that apparently dominates the matter density of the Universe. 
This is probably dominated by cold dark matter, with a density that probably falls within 
the range 0.2 < Q CDM < 0.5 [0, and may be in the form of massive weakly-interacting 
particles. It is therefore particularly important to search for such dark matter particles [||, 
and one of the most direct strategies is the search for relic particle scattering on nuclei in 
a laboratory detector |J. Many experiments around the world are engaged in this search, 
largely motivated by the cross sections calculated assuming that the cold dark matter is 
dominated by the lightest neutralino x @ °f the minimal supersymmetric extension of the 
Standard Model (MSSM) f§. 

We recently re-evaluated || the spin-dependent and spin-independent cross sections for 
neutralino scattering on protons and neutrons [0, assuming universality for all of the 
soft supersymmetry-breaking mass parameters of the MSSM including the Higgs multi- 
plets, incorporating the latest available LEP constraints on the MSSM parameter space, 
and assuming that the cosmological density of the relic neutralino falls within the range 
0.1 < Q x h 2 < 0.3, corresponding to the favoured range of flcDM and a Hubble expansion 
rate 0.6 < h < 0.8 in units Hq = 100 x h km/s/Mpc. We used the latest information from 
chiral symmetry || |10| , low-energy it —p,n scattering |ll| and deep- inelastic lepton-nucleon 
scattering |12] to fix the hadronic matrix elements. Our calculations fell considerably be- 
low the present experimental sensitivities [13|, as well as the highest theoretical estimates 
available in the literature ||14|| , some of which used less restrictive assumptions. There are, 
however, some other recent lower estimates: see |15[], for example, which is in good agreement 
with our previous work ||. 



Shortly after our paper appeared, the DAMA collaboration confirmed JT6| their previous 
evidence for the annual modulation of energy deposits in their scintillation detector, which 
they interpret as due to the scattering of some cold dark matter particle with mass between 
about 50 and 100 GeV, and spin-independent cross section on a proton between about 10~ 6 
and 10~ 5 pb. This cross section range is considerably larger than we found previously ||, 
though consistent with the range allowed by some previous cross section estimates. Sub- 
sequent to the DAMA paper, the CDMS collaboration has reported []T7| negative results 
from their experiment, establishing an upper limit on the spin- independent cross section 
that excludes most, but not all, of the range suggested by DAMA. 

This unresolved situation motivates us to explore more widely the possible neutralino- 
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proton cross sections in the MSSM, including both the spin-dependent and spin-independent 
(scalar) contributions. As before, we impose the latest constraints on the MSSM parameter 



space imposed by the LEP and other experiments [[18], |19|], such as measurements of b — > 57 
decay. It is important to note that the LEP limits we use here Jill have been updated 
significantly compared to what we used in ||. Notably, the chargino and particularly the 
Higgs mass limits we use here are stronger. The latter has a substantial effect at tan (3 = 3: in 
addition to the consequent direct reduction in the Higgs-exchange contribution to the scalar 
cross section, the improved lower limit on the Higgs mass further restricts mo and mi/2 from 
below, because of their contribution to rrih via radiative corrections. Also, previously we did 
not use the b — » 57 constraint, which we implement here by requiring rriA > 300 GeV for 
p < 0. All of these effects tend to remove some of the higher cross sections that we found 
previously, particularly at low m x . 

The main thrust of this paper, however, is to relax two of the theoretical assumptions 
made in our previous work. 

• The absence of large flavour- changing neutral interactions suggests that the soft super- 
symmetry-breaking scalar mass parameters m$ i of the MSSM may be universal for different 
quark and lepton flavours. However, there is no strong phenomenological or theoretical 
reason why the m Q . should be the same for the Higgs multiplets as for squarks and sleptons, 
and we relax this universality assumption in this work. It is known that, in this case, 
the lightest neutralino x might be mainly a Higgsino, but this particular option is greatly 
restricted by LEP data (IR M . 



• Neutralinos might not constitute all the cold dark matter, but might be complemented 
by other particles such as axions or superheavy relics. In this case, Q x < Qcdm, and 
Q x h 2 < 0.1 becomes a possibility. For any given neutralino mass, Q x may be decreased by 
increasing the x annihilation cross sections, which is often correlated with an enhanced elastic 
X-proton scattering cross section. Before concluding that cold dark matter detection becomes 
easier in this case, however, one must consider what fraction of our galactic halo density 
Phaio could be composed of neutralinos. Since the process of halo formation is essentially 
independent of the nature of the cold dark matter, as long as it is non-relativistic and weakly 
interacting, one should expect that 

P X = Phalo X (t^-) , (1) 

In an effort to be as optimistic as is reasonable, we assume that p x = phaio if ^\^ 2 — 0-1? 
and rescale: p x = phaio x (Q x h 2 /0.1) if Q x h 2 < 0.1. 

In our previous work 0, in which we assumed universality for the Higgs masses (UHM) 
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at the conventional supersymmetric GUT scale ~ 10 GeV, and the canonical range 0.1 < 
£l x h 2 < 0.3, we found that the possible ranges of elastic scattering cross sections were 
very narrow for any fixed values of m x ,tan/3 and the sign of /x, even allowing for plausible 
uncertainties in the hadronic inputs |Tt], [TT], |T2|] , and that they were always orders of 
magnitude below the present sensitivities |]13|| , even for the smallest allowed values of m x ~ 
50 GeV ||, H. Specifically, the maximum value we found for the spin-dependent x ~ P 



elastic scattering cross section for 3 < tan/3 < 10 was well below 10~ 3 pb, attained for 
m x ~ 60 GeV, and the maximum value we found for the spin-independent \ — p elastic 
scattering cross section for 3 < tan/3 < 10 was ~ 10~ 7 pb, again attained for m x ~ 60 GeV. 
The corresponding experimental sensitivities are ~ 1 pb and ~3x 10~ 6 pb, respectively. At 
higher neutralino masses, the predicted cross-sections were significantly smaller still. 

In the constrained version of the MSSM, when all soft scalar masses, including the Higgs 
masses, are set equal at the unification scale (UHM), there are four independent parame- 
ters, the soft scalar masses, mo, the gaugino masses, mi/2, the soft trilinear mass terms, A 
(assumed to be universal), and tan/3. In addition, there is the freedom to choose the sign of 
the Higgs mixing mass fi. Previously we scanned the mo — parameter space for fixed 
tan/3 and sgn(n). Our results were not very sensitive to A. 

Now that we relax the universal Higgs-mass assumption (nUHM) , we find much broader 
ranges of elastic scattering cross sections for any fixed values of m x ,tan/3 and the sign of 
/i. As previously, we perform a systematic scan of the region of the mo,mi/2 parameter 
space of the MSSM that is consistent with accelerator constraints. Here, mo refers only to a 
common squark and slepton mass, and the two Higgs soft masses m x and m 2 are fixed by the 
conditions of electroweak symmetry breaking, since we allow /i and the Higgs pseudoscalar 
mass rriA to be free parameters. Thus, we scan over mo, mi/2, /i, m^, and A for fixed tan/3. 
The details of these scans are given below, where we document which parameter choices fail 
which LEP constraint and/or the cosmological relic density requirement. 

We find that the elastic scattering cross sections may be somewhat larger than we found 
before in the UHM case, particularly for larger m x . However, the absolute values are still 
well below the present experimental sensitivities [[U|, at least for the canonical range 0.1 < 
f2 x /z 2 < 0.3 for the relic neutralino density. This remains true when we consider Q x h 2 < 0.1, 
but rescale the halo density as described above. 

We cannot exclude the possibility that there might be some variant of the MSSM that 
could accommodate the cold dark matter scattering interpretation of the DAMA data, but 
this would require an extension of the framework discussed here. One possibiliity might 
be to adopt a larger value of tan/3 [EDI: we restrict our attention to tan/3 < 10 to avoid 
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some uncertainties in the treatment of radiative corrections in the renormalization-group 
evolution of the MSSM parameters which affect the relic density calculations. Another 
possibility might be to relax further the universality assumptions for soft supersymmetry- 
breaking masses, either in the scalar or the gaugino sector. In particular, models in which 
rriq/mi is smaller than in the models discussed here might be able to accommodate larger 
elastic %-proton rates for any given value of Q x . Another way to reduce rriq/m^, with a 
similar effect, could be to postulate universality at a lower, intermediate renormalization 
scale, below the conventional supersymmetric GUT scale 

2 Theoretical and Phenomenological Background 

We review in this Section relevant aspects of the MSSM [||]. The neutralino LSP is the 
lowest-mass eigenstate combination of the Bino B, Wino W and Higgsinos Hip, whose mass 
matrix N is diagonalized by a matrix Z: diag(m Xll .. t 4) = Z*NZ~ 1 . The composition of the 
lightest neutralino may be written as 

X = Z xl B + Z x2 W + Z x3 Hi + Z xA H 2 (2) 

As previously, we neglect CP violation in this paper, so that there are no CP- violating phases 
in the neutralino mass matrix and mixing. For the effects of CP- violating phases on the neu- 
tralino scattering cross-section see f22fl-|25l. We assume universality at the supersymmetric 



GUT scale for the U{1) and SU{2) gaugino masses: M 1]2 = ^1/2, so that M x = | tan 2 9yyM 2 
at the electroweak scale. 

We also assume GUT-scale universality for the soft supersymmetry-breaking scalar masses 
m of the squarks and sleptons, but NOT for the Higgs bosons, in contrast to 0. We further 
assume GUT-scale universality for the soft supersymmetry-breaking trilinear terms A. Our 



treatment of the sfermion mass matrices M follows those in ]22| , and we refer the interested 
reader to pj for further details and notation. It suffices here to recall that, CP being con- 
served, the sfermion mass-squared matrix for each flavour / is diagonalized by a rotation 
through an angle 9f. We treat as free parameters mi/2 (we actually use M 2 which is equal to 
nii/ 2 at the unification scale), the soft supersymmetry-breaking scalar mass scale tuq (which 
in the present context refers only to the universal sfermion masses at the unification scale), 
A and tan (3. In addition, we treat fi and the pseudoscalar Higgs mass as independent 
parameters, and thus the two Higgs soft masses m 1 and m 2 , are specified by the electroweak 
vacuum conditions, which we calculate using m t = 175 GeV []. 

1 We have checked that varying m t by ±5 GeV has a negligible effect on our results. 
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The MSSM Lagrangian leads to the following low-energy effective four-fermion Lagrangian 
suitable for describing elastic x-nucleon scattering p2f : 

£ = Xl^XQilfiioiu + a 2 il 5 )qi + a 3 iXXQi1i + ^uxfx^Qi + a 5 iXXQa 5 Qi + otoiXI^XQiQi (3) 

This Lagrangian is to be summed over the quark generations, and the subscript i labels up- 
type quarks (i = 1) and down-type quarks {i = 2). The terms with coefficients 0^,0^,054 
and «6i make contributions to the elastic scattering cross section that are velocity-dependent, 
and may be neglected for our purposes. In fact, if the CP- violating phases are absent as 



assumed here, a 5 = a 6 = [23, 24]. The coefficients relevant for our discussion are: 
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where Di,T 3i denote hyper charge and isospin, and 



(4) 



(5) 



(6) 



hi — Z X 3(Z x4 ) , 5 2 i — Z x4 (—Z x3 ), 
Bi = sin/3(cos/9) , Aj = cos/5(— sin/?), 
Ci = sina(cosa) , Di = cos a(— sin a) 



(7) 
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for up (down) type quarks. We denote by m# 2 < the two scalar Higgs masses, and a 
denotes the Higgs mixing angle Q. 

As discussed in , the elastic cross section for scattering off a nucleus can be decomposed 
into a scalar (spin-independent) part obtained from the a 2 i term in (|3|), and a spin-dependent 
part obtained from the term. Each of these can be written in terms of the cross sections 
for elastic scattering for scattering off individual nucleons. We re-evaluated the relevant 
matrix elements in ||. Here we limit ourselves to recalling that: 

• There are uncertainties in the scalar part of the cross section associated with the ratios 
of the light-quark masses, which we take from |J: 

m m 

— = 0.553 ± 0.043, — = 18.9 ±0.8 (8) 

m d m d 

and information from chiral symmetry applied to baryons. Here the principal uncertainty 
is associated with the experimental value of the 7r-nucleon a term and the corresponding 
values of the ratios of the B q =< p\qq\p >. Following flQfl , we use 



Bu Bs = 1.49 (9) 
Bh — 



with a negligible experimental error, and [|TT 

2B 



B d + By 

which yields 



0.2 ±0.1, (10) 



— = 0.73 ±0.02 (11) 
B u 

The difference between the scalar parts of the cross sections for scattering off protons and 
neutrons are rather small. 

• The spin-dependent part of the elastic ^-nucleus cross section can be written in terms of 
axial-current matrix elements A^f^ that parametrize the quark spin content of the nucleon. 



We extract from a recent global analysis |12| the values 

Ai p) = 0.78 ± 0.04, Af = -0.48 ± 0.04, A® = -0.15 ± 0.04 (12) 

where the errors are essentially 100% correlated for the three quark flavours. In the case of 
the neutron, we have A&O = Ajf } , A^ n) = A<*\ and A^) = A<fl 



2 We note that (g, |) is taken from (23) and agree with @, g, || 
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3 Cosmological and Experimental Constraints 



As noted in ||, several convergent measures of cosmological parameters [IJ have suggested 
that the cold dark matter density VLcdm = 0.3 ± 0.1 and that the Hubble expansion rate 
H = h x 100 km/s/Mpc: h = 0.7 ± 0.1, leading to our preferred range 0.1 < QcDAih 2 < 
0.3. The recent data on the spectrum of cosmic microwave background fluctuations from 



BOOMERANG g6[ and MAXIMA |27j are consistent with this range, but do not signifi- 
cantly constrain it further. The upper limit on Qcdm can be translated directly into the 
corresponding upper limit on Q x . However, it is possible that there is more than one com- 
ponent in the cold dark matter, for example axions and/or superheavy relics as well as the 
LSP Xi opening up the possibility that Vt x < 0.1. For a given value of m x , values of the 
MSSM parameters which lead to Q x < 0.1 tend to have larger \ annihilation cross sections, 
and hence larger elastic scattering cross sections. Note, however, that the upper bound, 
Q x h 2 < 0.3, is a firm upper bound relying only on the lower limit to the age of the Universe, 
tu > 1-2 x 10 10 years (with f2 tota i < 1)- 

However, in such a 'shared' cold dark matter scenario, the packing fraction of neutralinos 
in the galactic halo must be reduced. As discussed in [P^j , for example, dark matter particles 
are taken into the halo in 'sheets' in phase space, whose thicknesses are determined by their 
initial (thermal) velocity. The 'sheets' of cold dark matter particles are of negligible thickness, 
so the ratios of their densities in the halo are identical with their cosmological densities, and 
therefore 

n O 

(13) 



Px _ ^x 



PCDM &CDM 

On the other hand, the 'sheets' of hot dark matter particles are of finite thickness related 
to their thermal velocities at the onset of structure formation, which limits the possible 
phase-space density of hot dark matter particles, so that Phdm / Pcdm < ^hdm/^cdm in 
general [p8fl . Moreover, a large ratio Qhdm I^cdm is currently not expected. 

The LSP detection rate also must be reduced correspondingly to (|i3|). Accordingly, 



when we consider MSSM parameter choices that have Q x h 2 < 0.1, we rescale the calculated 
scattering rate by a factor Q x h 2 /0.1. This rescaling by the minimal acceptable value of 
^CDAih 2 is relatively optimistic. 



For the calculation of the relic LSP density, we have included radiative corrections [18| to 



the neutralino mass matrix and include all possible annihilation channels ||29|| . In the MSSM, 
it is well known that there are large regions of the M2, /i parameter plane for which the LSP 
and the next lightest neutralino (NLSP) and/or chargino are nearly degenerate, namely in 
the Higgsino portion of the plane when M 2 ^> /i. It was shown |31| that, in these re- 
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gions, coannihilations between the LSP, NLSP, and charginos are of particular importance 
in determining the final relic density of LSPs, and these have been included in the present 
calculation. Inclusion of these coannihilation channels has the important consequence that, 
in the Higgsino regions where one expects larger elastic scattering cross sections, the relic 
abundance is substantially reduced. On the other hand, we do not include here coannihila- 
tions between the LSP and the sleptons t [32|, in particular the lighter stau fi, which were 
shown to play an important role in models with scalar mass universality also for the Higgs 
multiplets (UHM). These are known, in particular, to be important for determining the 
maximum possible generic value of m x in the UHM case, but are of less generic importance 
than x ~~ x' ~ coannihilations in the non-universal nUHM case considered here. For the 



same reason, we have also not implemented x ~ t coannihilation |j33[ . The neglect of such 
X — f coannihilation processes is generally conservative as far as the elastic scattering rates 
are concerned, since any reduction they cause in fl x h 2 is unlikely to be compensated by a 
corresponding enhancement in the elastic scattering cross section. We also do not pay any 
particular attention to the narrow parameter slice of mixed gaugino/Higgsino dark matter 
where oc mi/ 2 and m x may become large P3) , because this requires an adjustment of 
parameters at the % level, and is hence not generic. However, these are sampled, with the 
appropriate weighting, in our general randomized scan of the parameter space. 

The lower limit on m x depends on the sparticle search limits provided by LEP and other 
experiments |TR [B| . The most essential of these for our current purposes are those provided 



by the experimental lower limits on the lighter chargino mass m x ± and the lighter scalar Higgs 



mass m# 2 . As discussed in jyl, here we assume a lower limit m x ± > 101 GeV. The impact of 
the recently-improved lower limit on the Higgs mass ]35[ is potentially more significant [19 1 , 
particularly for tan/3 = 3, as displayed in Figs. 6 of fll~9l . The present experimental lower 
limit for tan/3 = 3 approaches m# 2 > 107 GeV [35). In implementing this constraint, we 



allow a safety margin of ~ 3 GeV in the MSSM calculations of m# 2 ||36|| , and hence require 
the MSSM calculation to yield m# 2 > 104 GeV for tan/3 = 3. In the case of tan/3 = 10, the 
LEP constraint on the MSSM Higgs mass is weaker (see Fig. 6 of [pjj), and we require only 
m# 2 > 86 GeV, which includes again a 3 GeV margin of uncertainty. The corresponding 



limit on m Q and in this case may be ignored fL9[ . The other two constraints that we 
implement are on sfermion masses, which we require to be (i) larger than 92 GeV, and (ii) 
larger than that of the lightest neutralino. We recall also the importance of the b — > 57 



constraint ||37|| , which we implement in an approximate way, by requiring > 300 GeV 
for jj, < |]l9fl . As also discussed in requiring our present electroweak vacuum to be 
stable against transitions to a lower-energy state in which electromagnetic charge and colour 
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are broken (CCB) [RSI would remove a large part of the cosmologically-favoured domain of 
MSSM parameter space. We have not implemented this optional requirement in the present 
study. In the next section, we will show the effect of the various expermental constraints on 
our scan of the parameter plane. 



4 MSSM Parameter Scan 

We have scanned systematically the MSSM parameter space, taking into account the cosmo- 
logical and experimental constraints enumerated in the previous Section and implementing 
the MSSM vacuum conditions for the representative choices tan/3 = 3 and 10. As discussed 
in [ 19 1 , lower values of tan /3 are almost entirely excluded by LEP. Our parameter scan was 






< m < 


1000 


80 


< < 


2000, 


80 


< M 2 < 


1000, 





< rriA < 


1000, 


-1000 


< A < 


1000. 



over the following ranges of parameters: 

(14) 
(15) 
(16) 
(17) 
(18) 

The main scan, which covers m ,/x and M 2 > 100 GeV and rriA > 300 GeV, was sup- 
plemented with smaller but significant subscans, to cover the smaller values of these four 
parameters as described below. The values of mo we use are fixed at the unification scale 
~ 10 16 GeV, while the values of the remaining parameters, //, il^m^, and A are evaluated 
at the electroweak scale. The lower cut off on both M 2 and /i is due to the lower limit on 
the chargino mass. As we indicated above, we impose a lower limit itla > 300 GeV for [i < 
to avoid problems with b — > 57. However, it should be noted that this restriction is quite 
conservative as, even for itla = 350 GeV, there are regions included in the above scan which 



are not allowed by b — ► 57 [19]. Similarly, even for \x > 0, where we impose no cut off on 
rriA, we have incuded some points which should be excluded on the basis of b — > sj. 

As can be seen in the Table, the overall scan was divided into three (four) specific regions 
for each value of tan/? and fi negative (positive), each with the number of points listed. 
The subscans with lower thresholds were designed to scour carefully the regions of MSSM 
parameter space close to the LEP exclusions, with the aim of ensuring that we sampled points 
close to their boundaries. For each subscan, we show the number of points which survive 
all the LEP experimental constraints discussed above, and we see that lower fractions of 
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Table V. Details of MSSM parameter scans, including the numbers of points that survive the 
LEP constraints and have a relic density in the favoured range. 



scan 


Total points 


survived LEP 


0.1 < tth 2 < 0.3 


M 2 ,n,m > 100, m A > 300 


tan (3 = 3, /i > 


30000 


17817 


1552 


tan (3 = 3, /i < 


30000 


17210 


901 


tan/5 = 10,/i > 


30000 


26498 


2588 


tan/3 = 10,/i < 


30000 


26507 


2337 


100 > M 2 ,/i > 80, m > 100, m A > 300 


tan f3 = 3, /i > 


20000 


75 





tan /3 = 3, /i < 


20000 


4410 


30 


tan/5 = 10,/i > 


20000 


1632 


14 


tan (3 = 10,/i < 


20000 


4480 


58 


M 2 , /i > 80, m < 100, m A > 300 


tan (3 — 3, /i > 


20000 


2669 


663 


tan /? = 3, /i < 


20000 


2247 


487 


tan/5 = 10,/i > 


20000 


5394 


2436 


tan/5 = 10,/i < 


20000 


5140 


2377 


M 2 , /i > 80, m > 0, m A < 300 


tan f3 = 3, /i > 


20000 


2208 


164 


tan/5 = 10,/i > 


20000 


12096 


1170 



the low-threshold subscans survive them, in particular because they tend to yield excluded 

values of the chargino mass. Fig. |] provides some insight into the impacts of the different 

LEP constraints for the case tan (3 = 3 and fi > 0. We plot in Fig. [I] the points scanned 

in the M 2 — fi parameter plane. In making this scatter plot, we show a randomly chosen 

subset of 5000 of the 90000 points sampled Q, since it is much easier to pick out the relevant 

physical effects of the cuts in such a subset of points, the full plot being extremely dense. 

We see that the chargino cut removes points at low values of /i and M 2 , denoted by 

(green) pluses, that the Higgs cut then removes many more points with low M 2 , denoted 

by (red) crosses, that the sfermion cut removes still more points with low M 2 , denoted 

by (violet) triangles (this occurs at high A and/or /i when there is a sizeable off-diagonal 

component in the sfermion mass matrix), and that the LSP cut tends to remove points at 

higher M 2 denoted by (golden) diamonds. The surviving (blue) squares are spread over the 

/i, M 2 plane, except for small values. Note that some points may fail to survive more than 

3 We have checked that there is no qualitative difference between this plot and the much denser plot with 
all points shown. 
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one of the above cuts. These are only denoted by the first cut tested and failed in the order 
listed above. The scans for the opposite sign of /i and for tan f3 = 10 exhibit similar features, 
and are omitted here. The only noticeable difference when \x > is that not so many were 
points eliminated by the Higgs cut at large values of M 2 for \i < 0, because we imposed 
the limit > 300 GeV: for fi > 0, many more points were run with low rriA and hence 
low rriH 2 - Also, for tan/3 = 10, more points survive at low \i and/or M 2 due to the relaxed 
contraint on the Higgs mass. 



tan {3 = 3, ^>0 




SO 10Q0 2000 



\x 

Figure 1: Results of the scan of MSSM parameter space for tan (3 = 3 and fi > summarized 
in Table 1, illustrating the impacts of the various LEP constraints. We denote by (green) 
pluses the points that fail the chargino cut: m x ± > 101 GeV, by (red) crosses the remaining 
points that survive the chargino cut but fail the Higgs cut: m H2 > 104 GeV, by (violet) 
triangles the points surviving the previous cuts that fail the sfermion cut: mj > 92 GeV, 
and by (golden) diamonds the points surviving the previous cuts that do not have the lightest 
neutralino as the LSP. The (blue) squares denote scan points that survive all these LEP cuts. 

The last column of the Table shows how many of the points that survive the LEP con- 
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Figure 2: Results of the scan of MSSM parameter space for tan (3 = 3 and // > summarized 
in Table 1, illustrating the impact of the cosmological relic density constraint on the points 
that points that survived the LEP constraints illustrated in Fig. [J. We denote by (green) 
pluses the points that have too small a relic density: Q x h 2 < 0.1, by (red) crosses the points 
that have too high a relic density: Q x h 2 > 0.3, and by (blue) squares the good points for 
which 0.1 < tt x h 2 < 0.3. 

straints and have relic densities in the cosmologically preferred range 0.1 < Q x h 2 < 0.3. It 
is apparent that most of the preferred points emerge from the first scan with M 2 , \i > 100, 
as the lower values which were explored thoroughly in the second subscan generally failed 
the chargino cut. More details of the scan over cosmological relic densities for tan/3 = 3 and 
fi > are shown in Fig. [|. As in Fig. |l], we show only a randomly selected subset of 5000 
points out of the total of approximately 22000 points which survived the LEP cuts. We see 
that, among the points that survived the previous LEP constraints, those with a small ratio 
of fx/M 2 generally have too small a relic density, denoted by (green) pluses, as a result of 
over-efficient x ~ x' ~ X coannihilation, whereas points with 11/M2 ~ 1 to 5 tend to have 
too large a relic density, denoted by (red) crosses, particularly if /j and M2 are individually 
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large. The points with a relic density in the preferred range 0.1 < Q x h 2 < 0.3, denoted by 
(blue) squares, tend to accumulate around fi/M^ ~ 1/2 or low M 2 . The former points are 
in the transition region between over-efficient \ ~ x' ~ X coannihilation and under-efficient 
annihilation at large fi and M 2 , whereas the latter are in the region of low M 2 where careful 
implementation of the LEP constraints is essential. However, it is apparent from Fig. |2] that 
there are exceptions to these general trends. We do not discuss them in detail, but remark 
that we have made an attempt to understand at least those exceptions that lead to 'unusual' 
elastic scattering cross sections. 

5 Elastic Scattering Cross Sections 

We now discuss the values of the elastic scattering cross sections that are attainable, bearing 
in mind the LEP and cosmological relic density constraints. Fig. |3] illustrates the allowed 
ranges of elastic scattering cross sections for the points included in our scan for the particular 
case tan (3 = 3, // > 0, as it was described in the previous Section. Plotted is a subset of 3000 
of the 90000 points scanned, indicating which points survive all the LEP cuts, and which 
other points fail which LEP cut. We find similar results for tan/? = 10 and/or the opposite 
sign of /i, with the exception that when p < we find some points trickling down below 
the apparent boundary at ~ 10~ 10 pb in Fig. 0(b), because of cancellations similar to those 
discussed in ||. 

We note, in particular, that the LEP chargino and Higgs cuts remove many points with 
low m x and/or large elastic scattering cross sections. The sfermion mass cut is less important. 
The constraint that x be the LSP removes quite a large number of points, populated more 
or less evenly in these cross section plots. The somewhat sparse set of points with very small 
cross sections give some measure of how low the cross section may fall in some special cases. 
These reflect instances where particular cancellations take place, examples of which were 
discussed in 0, and should not be regarded as generic. The lower boundary of the densely 
occupied region in Fig. ^ offers an answer to the question how low the elastic scattering cross 
sections may reasonably fall, roughly a ~ 10~ 9 pb for the spin-dependent cross section and 
~ 10~ 10 pb for the spin-independent cross section. 

We would like to draw particular attention to the spin-independent cross-section shown in 
Fig. 0(b). Notice that there are parameter choices with very large scattering cross sections. 
In this random selection, the cross-section may be as high as a few x 10~ 4 pb, and could 
even be larger than that claimed by DAMA. Indeed, in the full set of 90000 points scanned, 
there are even a few points which surpass 10 -3 pb. However, all of these points have been 
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Figure 3: Scatter plots of (a) the spin- dependent and (b) the spin-independent elastic scat- 
tering cross sections for tan /3 = 3, // > for a representative subsample of 3000 points, 
illustrating the impacts of the LEP constraints. The (green) plus signs denote points that 
fail the chargino mass constraint, which usually have small m x and sometimes large cross 
sections. The (red) crosses denote surviving points that fail the Higgs mass constraint, some 
of which have large spin-independent cross sections. The (violet) triangles denote surviving 
points that fail our (approximate) sfermion mass cut. The (golden) diamonds denote points 
where \ is not the LSP, and the (blue) squares denote points that survive all the LEP cuts. 

excluded by LEP (primarily by the Higgs mass cut). The largest surviving cross section is 
slightly over 10~ 7 pb, in both the randomly selected subset and the full scan. For fi < 0, 
the upper boundary in the scalar cross section is about an order of magnitude lower, as was 
the case in the model with universal Higgs masses ||. Note also that, for /i < 0, the limit 
rriA > 300 GeV we impose removes the points with large cross sections (in this case with 

^scalar ^ 10~ 8 pb). 

The next step is to implement the cosmological relic density constraints. We show in 
Fig. |] the cross sections obtained for a representative subsample of points with tan/3 = 
3,/i > that survive the LEP cuts, sorted according to the calculated values of Q x h 2 . 
Spin-dependent cross sections are plotted in panels (a) and (c), and spin-independent cross 
sections are plotted in panels (b) and (d). We include in panels (a) and (b) the cross 
sections calculated for unrealistic models with Q x h 2 > 0.3, and without making any rescaling 
correction for points with fl x h 2 < 0.1. The over-dense points with fl x h 2 > 0.3, denoted by 
(red) crosses, have been removed in panels (c) and (d), and the cross sections for under-dense 
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points with fl x h 2 < 0.1, denoted by (green) pluses, have been rescaled by the appropriate 
halo density fraction ([!]). As could be expected, the over-dense points tend to have smaller 
cross sections, and the under-dense points larger cross sections before applying the rescaling 
correction. After rescaling, the under-dense points yield cross sections in the range found 
for the favoured points with 0.1 < Q x h 2 < 0.3, denoted by (blue) boxes. For tan/5 = 10 
and n > 0, the scalar cross section is about an order of magnitude higher for points which 
survive all cuts. Relative to the cases with fi > 0, the < cases have a scalar cross section 
which is 1-2 orders of magnitude smaller. 

A comparison with Fig. |2| shows that the largest cross sections displayed in Figs. |](a,b), 
are almost all for Higgsino-like states whose elastic cross section is mediated by Z exchange. 
These are cosmologically under-dense, due to a combination of large annihilation and coan- 
nihilation cross sections. The cosmologically over-dense regions with relatively low elastic 
cross sections are mainly for gaugino-like states, and are for the most part more massive 



than 300 GeV, which is the oft-quoted upper bound on the bino mass in the MSSM [39]. 

Our resulting predictions for the spin-dependent elastic neutralino-proton cross section 
for tan (3 — 3 and \x > 0, after taking into account the LEP and cosmo logical constraints, are 
shown in Fig. |5](a), where a comparison with the UHM case is also made The raggedness 
of the upper and lower boundaries of the dark (blue) shaded allowed region reflect the 
coarseness of our parameter scan, and the relatively low density of parameter choices that 
yield cross sections close to these boundaries. We see that, at low m x close to the LEP 
limit, the spin-dependent cross section may be as much as an order of magnitude greater 
than in the UHM case considered previously ||, shown by the concave (red and turquoise) 
strip. However, even for low m x , the attainable range is far below the present experimental 
sensitivity, which is to a sp i n ~ 1 pb, and could be many orders of magnitude lower. As m x 
increases, the maximum allowed value of a spin decreases, though not as rapidly as in the 
previous UHM case |J. The hadronic uncertainties are basically negligible for this spin- 
dependent cross section, as seen from the light (yellow) shading. Turning now to the option 
tan (3 — 3 and \i < shown in Fig. ||(b), we see that the allowed range of the spin-dependent 
cross section is similar to that in the tan/3 = 3,/i > option. This is in contrast to the 
situation in the UHM || , where the spin-dependent cross section at low m x is much smaller 
for n < than for /i > 0. However, the cross section is still three or more orders of magnitude 
away from the present experimental upper limit. In the option tan/? = 10 and fi > shown 
in Fig. |5|(c), we see that the attainable range of the spin- dependent cross section is again 
similar to the previous option. This again contrasts with the UHM case, where the narrow 
4 In contrast to Q, here we have taken into account the updated LEP constraints. 
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allowed band for large m x ~ 500 GeV was somewhat higher than for the option tan/3 = 3 
and n > 0. As shown in Fig. ^|(d), our results for tan (3 = 10 and /i < are very similar to 
those for /i > 0. 

The analogous results for the spin-independent elastic neutralino-proton cross section, 
after taking into account the LEP and cosmological constraints, are shown in Fig. |6], where 
comparisons with the UHM case are also made. We see in panel Fig. ||](a) for tan/3 = 3 
and /i > a pattern that is similar to the spin-dependent case. For small m x , the spin- 
independent scalar cross section, shown by the dark (blue) shaded region, may be somewhat 
higher than in the UHM case, shown by the (red and turquoise) diagonal strip, whilst it could 
be much smaller. For large m x , the cross section may be rather larger than in the UHM case, 
but it is always far below the present sensitivity. The case shown in panel (b) of tan /3 = 3 
and fi < is somewhat different: the cross section never gets to be significantly larger than 
the UHM value at small m x . The reason for the anomalous extension of the UHM band 
outside the more general range is that the newer analysis reflected in the (blue and yellow) 
shaded region incorporates updated LEP constraints [HJ, that are significantly stronger for 



small tan/3 and small m x than those used in ||. This 'anomaly' is absent in panel (c) for 
tan/3 = 10 and /x < 0, which closely resembles panel (a), and also panel (d) for tan/3 = 10 
and /j > 0. We note in panel (d) a lesser reappearance of the 'anomalous' outdated UHM 
region at small m x . The dip in the (red and turquoise) UHM band for m x ~ 230 GeV in 
panel (d) reflects rather special cancellations |3| that are absent in the more general case. 
Overall, we note that the hadronic uncertainties, denoted by the light (yellow) bands, are 
somewhat larger in the spin-independent case than in the spin-dependent case. 

6 Summary and Prospects 

In this paper we have extended the analysis of || to consider a more general sampling of 
super symmetric models, relaxing the UHM assumption we made previously. For each of two 
choices of tan/3 and /i negative (positive), we have sampled 70000 (90000) sets of MSSM 
parameters, 30000 in general scans and 20000 each in two (three) special subscans over lower 
values of M2, fi, rriQ (and mi). We have implemented the current LEP constraints on MSSM 



parameters ||19| , discussing in detail which scan points survive which of these constraints. We 
have further discussed which of the remaining scan points yield a cosmological relic density 
in the allowed range Q x h 2 < 0.3, and which of these are in the preferred range Q x h 2 > 0.1. 
We exclude from further consideration the over-dense points with Q x h 2 > 0.3, and rescale 
the predicted cross sections for under-dense points with fl x h 2 < 0.1 as in ([!]). 
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The cross sections we predict for spin-dependent and spin-independent elastic neutralino- 
proton scattering for different values of tan (3 and the sign of fi are shown in Figs. |5] and [| 
respectively. We provide in Fig. [?| a compilation of our results, compared with the present ex- 
perimental upper limits on the cross sections [Tj| and the detection of spin-independent scat- 
tering reported by the DAMA Collaboration [IB|. We see that our predicted cross sections 
are well below the experimental upper limits for both the spin-dependent and -independent 
cases. We are unable to find MSSM parameter sets consistent with our relaxed universal- 
ity assumptions that come close to explaining the DAMA measurements. Our assumptions 
would need to be questioned if the neutralino scattering interpretation of the DAMA data 
is confirmed. For example, we have restricted our attention to models with tan/3 < 10. 
Alternatively, the DAMA data might favour models with values of irtq/m^ obtained either 
by relaxing the input universality assumption, or by imposing it at some renormalization 



scale below the conventional supersymmetric GUT scale [21 



In the future, we plan to improve the available relic density calculations by extending 
them to larger tan/3 and incorporating consistently all coannihilation processes. On the 
experimental side, we expect that other Collaborations will soon be able to confirm or 
exclude definitively the DAMA interpretation of their annual modulation signal as being 
due to neutralino scattering. Looking further ahead, we interpret our results as indicating 



a high priority for a new generation of direct dark matter detection experiments [ID with a 
much higher sensitivity. 
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Figure 4: Scatter plots (a,c) of the spin- dependent and (b,d) of the spin-independent elastic 
scattering cross sections for tan/3 = 3, fi > 0, after implementing the LEP constraints, 
exhibiting the impacts of the cosmological relic density constraints. The (green) pluses have 
Q x h 2 < 0.1, the favoured (blue) boxes have 0.1 < Q x h 2 < 0.3, and the (red) crosses have 
Q x h 2 > 0.3. Note in panels (c) and (d) the impacts of removing the over-dense points, which 
tend to have lower cross sections, and rescaling the under-dense points as in ([|), suppressing 
some high cross section points. 
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Figure 5: Allowed ranges of the spin- dependent elastic neutralino-proton cross section for 
(a) tan/3 = 3 and fi > 0, (b) tan/3 = 3 and fi < 0, (c) tan/3 = 10 and fi > and (d) 
tan/3 = 10 and fi < 0. The main (blue) shaded regions summarize the envelopes of possible 
values found in our scan, for points respecting the LEP constraints, discarding points with 
fl x h 2 > 0.3, and rescaling points with Q x h 2 < 0.1 according to (0). The small light (yellow) 
shaded extensions of this region reflect the hadronic matrix element uncertainties discussed 
in Section 2. The concave (red and turquoise) strips are those found previously assuming 
universal Higgs scalar masses (UHM) 
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Figure 6: Allowed ranges of the spin-independent elastic neutralino-proton cross section for 
(a) tan/3 = 3 and fi > 0, (b) tan/3 = 3 and \i < 0, (c) tan/3 = 10 and fi > and (d) 
tan/3 = 10 and fi < 0. The main (blue) shaded regions summarize the envelopes of possible 
values found in our scan, for points respecting the LEP constraints, discarding points with 
fl x h 2 > 0.3, and rescaling points with £l x h 2 < 0.1 according to (0). The small light (yellow) 
shaded extensions of this region reflect the hadronic matrix element uncertainties discussed in 
Section 2. The (red and turquoise) diagonal strips are the results found assuming universal 
Higgs scalar masses (UHM) j^J. 
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Figure 7: Compilations of our allowed ranges for (a) the spin- dependent elastic neutralino- 
proton cross section, and (b) the spin-independent elastic neutralino-proton cross section for 
both the values of taxi /3 and the signs of fi studied. Our results are compared in panel (a) 
with the available experimental upper limits [15j, and in panel (b) with the detection reported 
by the DAMA Collaboration J76|/ ; as well as with upper limits from other experiments [71 
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